Quasiperiodic Heisenberg antiferromagnets in two dimensions 
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We describe some of the properties of 2d quantum quasiperiodic antiferromagnets as reported in 
studies that have been carried out in the last decade. Many results have been obtained for perfectly 
ordered as well as for disordered two dimensional bipartite quasiperiodic tilings. The theoretical 
methods used include spin wave theory, and renormalization group along with Quantum Monte Carlo 
simulations. These methods all show that the ground state of these unfrustrated antiferromagnets 
have Neel type order but with a highly complex spatial distribution of local staggered magnetization. 
The ground state properties, excitation energies and spatial dependence, structure factor, and local 
susceptibilities are presented and discussed. The effects of introducing geometrical disorder on the 
magnetic properties are discussed. 

PACS numbers: 75.10.Jm, 71.23.Ft, 71.27.+a, 75.10.-b 



I. INTRODUCTION 



This article reviews properties of the Heisenberg model 
for a quantum spins in a quasiperiodic antiferromagnet. 
The quasiperiodic structures considered are two well- 
known models of quasicrystals, the Penrose and the 
Ammann Q tilings. Quasicrystals are highly organized 
states of matter (see with long range positional order 
that is a generalization of the usual translation invariance 
found in crystals. In a crystal, one can identify a basic 
motif which is replicated at regular distances, whereas in 
a quasicrystal, one can discern patterns which repeat in 
space, in aperiodic fashion. The simplest example of an 
aperiodic structure is a one-dimensional array of atoms 
with two characteristic spacings, say L and S (for long 
and short respectively) in which the spacings follow the 
Fibonacci sequence, ...LSLSLLSLS.. .according to a per- 
fectly deterministic algorithm. Due to the complexity of 
aperiodic structures compared with crystals, they have 
sometimes been considered as being in some sense closer 
to disordered structures insofar as their physical proper- 
ties are concerned. As we will see, on the contrary, the 
symmetries present in quasicrystals can lead to novel and 
striking features that are not found in a crystal, much less 
a disordered material. 

The presence of Bragg peaks in the structure fac- 
tor is one defining characteristic of a perfectly ordered 
quasiperiodic structure. Unlike crystals, where Bragg 
peaks occur at regular intervals determined by a set of 
d reciprocal lattice vectors (where d is the dimension- 
ality), the Bragg peaks of a d-dimensional quasicrystal 
are described with the help of n > d reciprocal lattice 
vectors. The peaks are not regularly spaced, and their 
intensities cannot be simply expressed as a function of 
the position. In principle, the reciprocal space is densely 
filled with peaks most of them of vanishingly small in- 
tensity, so that in an actual experiment, one sees a large 
but not infinite number of peaks that have an amplitude 
greater than some experimentally observable threshold 



value. For real quasicrystal samples, moreover, the ex- 
perimentally observed diffraction spots have finite widths 
due to instrumental resolution and sample imperfections. 
The widths of the spots indicate the correlation length, 
or distance out to which quasiperiodic order is present, 
and this is comparable to those observed in crystalline 
materials, of the order of microns in the best samples 
available. Since the constraints imposed by translational 
invariance are not present for the quasicrystal, the usual 
crystallographic rules forbidding certain rotational sym- 
metries such as 5-fold symmetry, do not apply. The ob- 
servation of five-fold symmetry of the diffraction pattern 
led to the experimental discovery by Schechtman et al Q 
of the first quasicrystalline alloy, composed of aluminum, 
copper and silicon, in 1984. Other crystallographically 
forbidden symmetries, such as 8-fold, 10- fold or twelve- 
fold symmetries have also been discovered and known 
quasicrystal alloys now number in the hundreds. 

In the above we have assumed that the structure is per- 
fectly ordered, however, many quasicrystals are thought 
to be intrinsically disordered, and better described by 
random tiling models rather than the above-mentioned 
deterministic structures. Random tilings [6j are com- 
posed of the same units as the deterministic tilings, but 
do not have the hierarchical invariance of the perfect sys- 
tems. A perfect structure can be disordered, for example, 
by allowing for random local phason flips - when an atom 
has the freedom to flip between two alternative positions 
without affecting the local bond lengths and angles. Here 
we are interested in the effect of this type of disorder on 
the antiferromagnetic ground state properties. 



A. Magnetic quasicrystals and the Heisenberg 
model 



The Heisenberg model can be used to discuss quasiperi- 
odic alloys such as the ZnMgR alloys (where R stands for 
a rare earth Ho,Ti,Dy,..), where the magnetism is due to 
localized f-orbitals. Properties of transition metal alloy 
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systems such as AlMnPd which have itinerant magnetic 
moments will not be discussed here (see review by Hip- 
pert et al in Q)- Magnetization measurements for the 
icosahedral ZnMgR and CdMgR compounds show that 
the predominant interactions are antiferromagnetic 
11 1 . Neutron scattering experiments 12-1^ show that 
although short magnetic correlations set in at low tem- 
perature, there is no long range antiferromagnetic or- 
der down to the lowest temperature studied. Instead, 
in some of these compounds, a spin-glass type low tem- 
perature phase and out-of-equilibrium phenomena was 
found. These phenomena arise due to magnetic frustra- 
tion as well as, no doubt, to the presence of some amount 
of disorder. Crystal field effects have been argued to be 
important as well [l^. The experimental data raise a 
number of interesting theoretical questions. What is the 
ground state of a quasiperiodic antiferromagnet and how 
is it different from that of periodic antiferromagnets ? 
What is the equivalent of spin wave modes or magnons 
in a quasicrystal ? What are the effects of disorder ? Are 
there novel frustration-induced phenomena ? 

These questions are difficult to answer, especially as 
regards the effects of disorder and frustration, whose ef- 
fects remain insufficiently understood even in periodic 
structures, despite the many studies devoted to them. 
We will here only consider the simplest situation: that 
of systems without frustration - those in which the spins 
have a ground state with a Neel type (up-down) order- 
ing. The most tractable of quasiperiodic models is the 
Fibonacci chain, and the behavior of quantum spins on 
quasiperiodic chains been considered by several authors. 
Quantum spin chains have been analyzed using renormal- 
ization schemes [11-13 based on the inflation symmetry 
of these chains. Using a mapping to fermionic models and 
techniques of bosonization, results have been obtained for 
global properties such as the magnetization as a function 
of external field, and the spectral gaps for a variety of 
different quasiperiodic sequences [19|] . In contrast to two 
and higher dimensions, however, there is no symmetry 
breaking in one dimension. 

Two and three dimensions are interesting because long 
range magnetic ordering can occur at zero temperature 
(destroyed at finite temperature in the case of d=2). 
In two dimensions, the low connectivity leads to larger 
quantum effects, compared to three dimensions. In addi- 
tion, some analytical calculations are possible for two di- 
mensions, and numerical calculations can be more easily 
carried for large systems compared to three dimensional 
structures. We will be interested in the quantum fluctu- 
ations in the quasicrystal, which are maximal when the 
spin S of the atoms is small. We will first describe the 
antiferromagnetic ground state of perfect quasicrystals in 
some detail, before going on to the case of disordered qua- 
sicrystals. It should be noted that this quantum limit is 
not expected to be pertinent in the rare earth quasicrys- 
tals where the spin S is on the contrary quite large, and 



the classical Heisenberg model should suffice. Competing 
nearest-neighbor and next-nearest neighbor interactions 
in these systems would lead in these systems to complex 
behavior at low temperatures. Such models are outside 
the scope of the current review, but have been studied 
for classical spin systems in two dimensions (20l-24 and 



three dimensions [15|, |25 1 . 



In recent years, new progress has been made in the 
attempt to fabricate quasicrystalline structures by de- 
position of atoms on clean surfaces of perfect quasicrys- 
tals. This is an exciting development, as two dimen- 
sional structures based on one or two adatom species 
are easier to model theoretically, making it possible to 
study in a controlled way the consequences of long range 
quasiperiodic geometric ordering on the electronic and 
magnetic properties. There has been much recent work 
on growing thin films on quasicrystalline surfaces that 
display quasiperiodic order out to distances of hundreds 
of nanometers. Characterization of these structures, and 
studies of the interplay of geometrical and electronic 
properties of these films are ongoing ( as reviewed in [26|). 
Magnetic properties of atoms in surface layers have not 
yet been investigated, and the realization of a Heisenberg 
model in thin quasiperiodic films represents an experi- 
mental challenge. 

Sec. II presents an overview of the two tilings with some 
definitions and geometrical properties. Sec. HI provides 
an introduction to properties of the Heisenberg model 
for quantum spins. Sec. IV introduces the models that 
are studied, whose results are then discussed in Secs.V 
and VI. Sec. VII discusses disordered models. Sec. VIII 
presents the conclusions. 



II. TILING MODELS 

Investigations of quasicrystal properties often begin 
with studies based on quasiperiodic tilings. In the way 
that the Bravais lattices are used as templates for crys- 
tals, quasiperiodic tilings are useful as a starting point 
for a description of real alloys. Depending on the dimen- 
sionality and symmetries of the tilings, different tilings 
help to classify the various members of this large class 
of materials. In three dimensions, the icosahedral tiling 
(also called 3D Penrose tiling) [27H30[ is the basis for 
most models of realistic atomic structures of quasicrys- 
tals. In one dimension, the simplest example of a tiling is 
the Fibonacci "tiling" or chain already referred to above. 
The best known tiling is, probably, the 2D Penrose tiling 
[H 0] , which was introduced before the discovery of qua- 
sicrystals. It is built out of elementary rhombus-shaped 
tiles that come in two shapes - the thick and the thin - 
each of which can occur in five orientations as can be seen 
in Fig. 1, leading to the five- fold symmetry of the diffrac- 
tion pattern. This structure gives rise to delta-function 
peaks in its structure factor (the modulus squared of the 



Fourier transform of a crystal formed by, say, placing a 
point mass on every vertex) . The octagonal or Ammann- 
Beenker tiling is another often studied two-dimensional 
tiling. It is composed of two kinds of tiles: the square 
and the 45° rhombus - the tiles shown in Figl^l 




FIG. 1: A finite patch of the perfect Penrose tiling. 




FIG. 2: Inflation of the octagonal tiling showing the original 
tiling (thin lines) and the inflated one (thick lines) 



A. Symmetries 

One readily observable feature of a tiling such as the 
one in FiglT]is that, despite the absence of perfect trans- 
lational invariance, patterns of any given size frequently 
repeat (this "local isomorphism" property is explained by 
Baake in (sij). The mean repetition distance of a pat- 
tern of linear size R is of the order of R, a property that 
replaces the strict translational invariance of crystalline 
structures. 

Next, we discuss the rotational symmetry present in 
these tilings. This holds in the same weaker sense than 
in crystals, as we saw with the repetitivity property. The 
rotational invariance is true in the following sense - any 
given configuration of tiles of arbitrary size R will be 
found to occur with equal probability in each of its n 
rotated forms, where n = 5 for the Penrose rhombus 
tiling and rt = 8 for the octagonal tiling. 

The tilings possess a hierarchical symmetry under dis- 
crete scale transformations. These are the so-called in- 
flation and deflation transformations illustrated in Figl2l 
Inflation (deflation) is a reversible operation of change of 
scale which consists of combining (subdividing) each of 
the tiles into larger (smaller) tiles of the same shape. In 
the case of the Penrose tiling, the tiles of the new system 
thus obtained are bigger (smaller) by a factor r = 
(the golden mean). In the case of the octagonal tiling, 
the edge length is increased by the number A = -I- 1 
(which goes by the name of silver mean). No new envi- 
ronments are created or destroyed in the process of in- 
flation or deflation. The original and the inflated tilings 
are equivalent in the sense that for any given configura- 
tion of tiles in one of the tilings, one can always find an 
exactly identical region (up to a length change) in the 
other tiling. 

The self-similarity of quasiperiodic tilings under 
inflation-deflation operations is an important symmetry 
property, and can be used in defining renormalization 



transformations. In the one-dimensional case, this prop- 
erty has been exploited to calculate properties of spin 
models on the Fibonacci chain and other so-called metal- 



lic mean chains [17[. There are fewer results available in 
the case of two dimensional tilings. A renormalization 
argument was used in [2^ to map out a phase diagram 
for a frustrated Ising model on the Penrose tiling. In the 
case of the octagonal tiling, an approximate renormal- 
ization group calculation leads to a nontrivial fixed point 
solution for the antiferromagnet ic g round state, as will 
be discussed in the next section [32, l33 1 . 

The last symmetry property that we will mention has 
to do with the bipartite character of the tilings. Quadri- 
lateral (or more precisely, rhombus) based structures 
such as the ones shown in FigH] are constituted of two 
sublattices. Less obviously, it can be shown that for the 
infinite quasiperiodic tiling, the two sublattices are ex- 
actly equivalent, as regards the number and the nature 
of the vertices (see below). 



B. Local environments 

The upper figure in FiglS] shows the different kinds 
of vertices or local environments present in the Penrose 
rhombus tiling. The seven local environments correspond 
to site coordination numbers (number of nearest neigh- 
bors) ranging from 3 to 7. The frequencies of each of 
the vertices can be calculated in terms of r. The average 
value J being exactly 4, as it is in the octagonal tiling. 
The lower figure shows a portion of the octagonal tiling 
labelling six different types of vertices, with the coordina- 
tion numbers ranging from 3 to 8. Again, the frequencies 
of each kind of vertex can be exactly found, in terms 
of the number A. The transformation of vertices under 
inflation or deflation on the Penrose and the octagonal 
tiling will not be described here in detail. We will return 
to this point in the section on renormalization group for 
the octagonal tiling. 





FIG. 3: (top) Local environments in the Penrose tiling (Re- 
produced from [sH])- (bottom)A portion of the octagonal 
tiling showing the six different nearest neighbor environments 
A,B,...,F (Reproduced from [33j) 



C. Tilings as projections from higher dimensions 

There are many ways to obtain quasiperiodic tilings. It 
is outside the scope of this review to describe in detail the 
methods that exist. In this section we will briefly outline 
the cut- and- project method used to obtain samples of the 
Penrose tiling in our calculations, since it will turn out 
to be helpful for understanding the results that will be 
obtained for the Heisenberg spin model. 

The Penrose tiling (PT) and octagonal tiling (OT) can 
be obtained by projecting a subset of vertices of a D di- 
mensional cubic lattice [3J] , where I? = 5 for the Penrose 
tiling and D = 4m the case of the octagonal tiling. In the 
projection method, the d-dimensional tiling is obtained 
by projecting a subset of points of a D-dimensional pe- 
riodic lattice (where D > d). The value of D and the 
selection of the subset of points to be projected are found 
by considering the symmetries one would like the tiling 
to possess. The Fibonacci chain, for example, can be 
obtained by projecting all the points lying within an in- 
finite strip of a square lattice of slope l/r, as shown in 
FigHl Inside the strip, whose width corresponds to one 
unit square, neighboring vertices are connected by edges 
whose projections give rise to the L and S line segments 
that make up the tiling. The figure shows by a thick 
line the set of edges that are selected for projection. The 
irrational slope of the strip ensures that the projected 
pattern never repeats. 

If, now, the irrational number were to be replaced 
by a nearby rational, for example one of the approxi- 
mants of the golden mean, 1/rfc = Fk/Fk+i, the strip 
becomes commensurate with the square lattice. It is easy 



FIG. 4: Obtaining a piece of the Fibonacci chain by the cut 
and project method. The bold line indicates the edges which 
are selected and form the Fibonacci chain after projection 
onto the physical line Ei. The perpendicular space axis is 
denoted £2. 



to check, using the recursion relation for the Fibonacci 
numbers, that the projected pattern would in that case 
repeat after a distance Lk = Fk+2 ■ This type of so-called 
approximant of the Fibonacci chain is useful for numer- 
ical calculations, where one considers finite systems with 
periodic boundary conditions, and attempts to extrapo- 
late to the quasiperiodic limit where fc — > 00. The infinite 
chain can be transformed into itself under inflation (or 
deflation) ~ in an inflation, tiles are grouped into new 
bigger tiles L' and S' {LS -> L' and L -> S"). The new 
chain obtained in this way has the same properties as the 
old, with the difference being a scale factor r. 

In an analogous way, one can obtain the octagonal 
tiling by projecting lattice points lying inside an infi- 
nite cylinder in four dimensional space. The orientation 
of the cylinder is given in terms of the so-called "silver 
mean" A = {V2 + 1). It is this irrational number that 
determines the shapes of the tiles after projection, and 
it determines the scale of the infiation/defiation trans- 
formation of the octagonal tiling. The Penrose tiling in 
two or three dimensions, similarly, may be obtained by 
projection, and the irrational number r appears again in 
the matrices used for projection and in the tile-rescaling 
transformations. As in the case of the Fibonacci chain, 
it is possible to obtain approximants of the octagonal 
tiling. These are finite square tilings that can be period- 
ically continued in the plane, and whose period can be 
increased progressively by considering successive approx- 
imants of A [35|. Similarly, for numerical studies of the 
Penrose tiling antiferromagnet, one can work with the so- 
called Taylor approximants obtained by a generalization 
of the method described by Duneau and Audier in [s^] • 

When the vertices of the tiling are projected not onto 
the physical plane, but instead into perpendicular space, 
they fall into a compact region (the "selection window"). 
In the case of the Fibonacci chain, one sees that points 
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project onto the line segment labelled by W in FiglH In 
the case of the octagonal tiling, all the projected points 
fall into an octagonal shaped domain in the two dimen- 
sional perpendicular space. The vertices of the Penrose 
tiling can be similarly projected into a three dimensional 
perpendicular space, where they lie in pentagon-shaped 
plane domains as seen in Fig ITT] In each of these cases, 
sites that have the same local environment are found 
grouped together in compact domains. In other words, 
the perpendicular space representation allows us to sep- 
arate sites not only according to their nearest neighbor 
configurations (coordination number) but also out to fur- 
ther distances. We will later exploit this useful property 
of the tilings in order to represent the complex spatial 
structure of the antiferromagnetic ground state in a com- 
pact way. 

The codimension, or the difference Z? — rf, is a rough 
measure of the degree of complexity of the quasiperiodic 
structure obtained in the projection method. The octag- 
onal tiling has a co-dimension of 2. The Penrose tiling 
can be obtained from a projection down from D = b (one 
can bring this down to Z? = 4, but with a less transparent 
procedure), and has a slightly larger set of local environ- 
ments than the octagonal tiling. If one projects a three 
dimensional cubic lattice onto a suitable plane one ob- 
tains the two-dimensional generalized Rauzy tiling [36| . 
which is of co-dimension 1. This last example is close 
to being an integrable system, as shown by the analyses 
of its spectral properties for the tight-binding electron 
problem [s^]- Going to higher dimensions and codimen- 
sions > 2, one has the entire family of rhombus tilings 
described in 38[. These have increasingly complex envi- 



ronments as the codimension gets larger. 



III. INTRODUCTION TO QUANTUM 
ANTIFERROMAGNETS 



The Heisenberg model expresses the energy as a func- 
tion of the orientation of the local moments of the atoms, 
and can be derived as an effective low energy model for 
interacting electron systems. The model considers a set 
of spins Si, situated on the sites z of a lattice or tiling and 
interacting via the Hamiltonian H — JijSi.Sj. The 
signs and the magnitudes of the spin-spin couplings Jij 
depend on the physical and chemical details of the struc- 
ture, and these variations give rise to the large variety of 
magnetic phases observed in solids. 

The simplest case, ferromagnetism, corresponds to one 
in which all of the couplings have negative sign Jij < 
- it is then energetically favorable to have spins aligned 
parallel to each other. This results in a unique ground 
state configuration (upto global rotations) . The resulting 
ground state has a simple description in classical as well 
as in quantum theory. This is not true in the antiferro- 
magnetic case, when the coupling J has a positive sign. 



as we now explain. To begin with, a new phenomenon 
called frustration can now occur. To illustrate this, let us 
consider antiferromagnets in which the coupling between 
nearest neighbor spins is J > 0. For simplicity, let us 
assume other couplings can be neglected. Unfrustrated 
systems are ones which comprise a sublattice A in which 
all spins point "up" (with respect to the axis of sym- 
metry breaking) and a sublattice B in which spins point 
" down" , in such a way that all neighbors of A-sublattice 
sites belong in the B-sublattice, and vice versa. FigE] 
shows one example of an unfrustrated and an frustrated 
loop. The unfrustrated loop has one possible minimal 
energy configuration (upto global transformations). The 
frustrated triangular loop minimizes its energy when the 
three spins are aligned at 120° angles with respect to 
each other, and this is possible in two different ways. In 
a large system formed from many such units, as a result, 
frustration can sometimes lead to very highly degenerate 
ground states, and such systems do not show a transition 
down to very low temperatures (T << J). 

Unfrustrated antiferromagnets are expected to have a 
unique ground state. For classical systems the ground 
state is characterized by broken rotational symmetry, 
with the spins ordering in an up-down "Neel" order. The 
order parameter, called the "staggered magnetization" , is 
defined by ms{T) = \Siz\ where z is the axis of symmetry 
breaking. In three dimensions, is zero above a critical 
temperature Tc, while below Tc, it grows as the temper- 
ature is decreased towards the maximum value possible. 
In one and two dimensions, ms{T) = at any nonzero 
temperature (Mermin-Wagner-Hohenberg theorem). 

At T = and for classical spins in any dimension, 
the order parameter is mf{T — 0) — S, where S is the 
spin length. In contrast, in the case of quantum spins, 
the order parameter rnf{T = 0) is reduced below this 
value, due to quantum fluctuations, which exist even at 
T — Q. Quantum fluctuations become more important 
as the spin length decreases, and are most important 
ioY S — i. The magnitude of the order parameter de- 
pends on the type of lattice, and it can vanish if the 
connectivity is sufficiently low, as happens in the one di- 
mensional chain [s^l- This is of course also true for the 
one-dimensional quasiperiodic systems: such as the Fi- 
bonacci chain Heisenberg antiferromagnet 3) i3]- Two 
dimensional models were an open question for a while. 
Interest in the two dimensional Heisenberg antiferromag- 
net rose to a peak in the mid-eighties with the discovery 
of the high-Tc superconductors^Many different analyti- 
cal and numerical approaches [40] have been brought to 
bear on the square lattice and other periodic antiferro- 
magnets. The challenges faced are due to the fact that 
unlike ferromagnets, the ground state of the antiferro- 
magnet is impossible to write explicitly for even a simple 
case such as the square lattice. There is close compe- 
tition for ground state status between the antiferromag- 
netic state and a state with no magnetic long range order 
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such as Anderson's RVB - for resonating valence bond - 
state 41| which resurfaced in the context of high- Tc su- 
perconductor physics (4^ . 

The results of these studies have settled the contro- 
versy for most of the community, and it is accepted that 
there is indeed long range antiferromagnetic order for 
spins 5 = ^ at T = in the two important unfrustrated 
two dimensional periodic lattices, namely the square and 
honeycomb. As coordination number decreases, quan- 
tum fluctuations increase. This is illustrated by the fact 
that the order parameter is smaller in the honeycomb 
lattice (tos = 0.235 [ill) compared to the square lattice 
(nis = 0.3070). 

Our aim here is to study the interplay of magnetism 
with the lattice geometry in an unfrustrated quasiperi- 
odic antiferromagnets in two dimensions by consider- 
ing the Heisenberg model on the octagonal and Pen- 
rose tilings. We will consider both perfect and phason- 
disordered (see last section) samples of these tilings. 
Other tilings in two dimensions that it would be interest- 
ing to study are the 2d Rauzy tiling 36] , which is closer 
to the square lattice than the previous two, and, at the 
opposite extreme in terms of complexity of structure, the 
high codimension random rhombus tilings introduced by 
Destainville et al 38|. The labyrinth [45'| tiling, which 
has the connectivity of the square lattice, and shares the 
same ground state properties if the couplings are assumed 
to be uniform in value. 



(I) 




(U) 



\ 



FIG. 5: Examples of minimal energy spin configurations: 
in an unfrustrated (left) square where the spins alternate in 
direction, and a frustrated (right) triangle where the energy 
is minimized by two different configurations I and II (A and 
B are two vectors oriented at ±120° angles with respect to 
the third spin). 



spatial complexity of the solutions of this model and non- 
trivial magnetic ground state and spin excitations arise, 
thus, solely due to the complexity of the quasicrystalline 
structure in real space. The Hamiltonian takes the form 



(id) 



I J Si- Sj 



(1) 



where Si is a spin variable on the site i, situated at po- 
sition ri. The antiferromagnetic couplings Jij = J > 
for all pairs of sites i and j that are linked by an edge, 
and zero otherwise. We will consider the case where the 
quantum number S — ■^j the case of greatest theoretical 
interest, as quantum effects are expected to be strongest. 
The systems we consider are bipartite, i.e. each site be- 
longs to one of two sublattices A and B, and the J-terms 
couple a spin on sublattice A to spins lying on sublattice 
B and vice- versa. This property ensures that there is no 
frustration, i.e., if one considers classical spin variables, 
the ground state is one for which all bonds are "satisfied" 
- with all the A-sublattice spins pointing in one direction 
and B-sublattice spins pointing in the opposite direction. 

This quasiperiodic unfrustrated Heisenberg antiferro- 
magnet is expected to have long range order at zero tem- 
perature, like the unfrustrated periodic two dimensional 
systems already mentioned. In the ground state, rota- 
tional symmetry is broken, and spins acquire a non-zero 
magnetization in the z-direction. In the quasicrystal, one 
can define a local staggered magnetization for every site 



(2) 



where Ci — ±1 depending on whether the site i is on the 
A or B sublattice. In a translationally invariant system, 
the local staggered magnetization is independent of i, 
while in the quasicrystal each site is different from all 
the others and we expect, in general, a distribution of 
values of nisi- Due to the symmetries of the quasicrystal, 
we will see that this distribution possesses very unique 
spatial properties. The methods that have been used to 
analyse the ground state properties are described below. 



A. The Heisenberg star cluster 



IV. MODELS OF QUASIPERIODIC 
ANTIFERROMAGNETS AND METHODS OF 
SOLUTION. 

We consider that spins are coupled only to their near- 
est neighbors, as defined by the edges of the tiling. We 
will assume all couplings to have the same value J. The 



We begin with a description of the energy of a simple 
cluster of spins, the Heisenberg stars. A central spin 
So is coupled to z neighbors Sj as shown in FiglHl The 
corresponding Hamiltonian, Hhs{z) = JJ2j=i^o ' ^jj 
can be exactly solved. The result for the lowest energy 
level of the cluster is 

e("\z) = ^JS{zS+l)^-/-^ (3) 
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where the second equahty holds for the case oi S = ^. 
The superscript denotes that this is the energy of an 
isolated cluster. This expression of the cluster energy 
is the first step, as we will explain later, for a recursive 
calculation of the energy in an extended structure. The 
degeneracy of the ground state manifold is 2Stot + 1 = 2;, 
rotational symmetry is not broken, and the value of the 
staggered magnetization is zero for all the spins. One can 
use the average value of the spin-spin correlation function 
to define a quantity that has the same dimensions as the 
magnetization squared, namely 



state energy. The value of the GSE found for the octag- 



^((So.So + ^|So-S,)o|) 
l-(5(5 + l)-6(")(z)/J) 



z + 1 



5 + z 
4(^ + 1) 



(4) 



where (..)o are the expectation values evaluated in the 
ensemble of ground states, and where the last equality 
holds for the case of 5* = ^. Note that for z large, /i 
tends to the classical value of the magnetization, fj.(z — > 
00) — >■ 5 and also that the quantum corrections to fi 
decrease with increasing z. 




FIG. 6: Heisenberg star with z neighboring spins. 



B. Quantum Monte Carlo calculations 

Relatively large - of the order of several thousand spins 
- samples can be studied using Stochastic Series Expan- 
sion Quantum Monte Carlo algorithm H CSSE-QMC). 
The first calculations were carried out [46| for the oc- 
tagonal tiling antiferromagnet. The ground state energy 
per site (GSE) of the quasicrystal was determined by re- 
peating calculations for a series of samples of finite size 
N (number of sites). The finite size scaling argument 
used for periodic systems predicts that the corrections 
to the ground state energy decrease as N~^^'^ for two- 
dimensional lattices [13] . The argument used can be gen- 
eralized to the quasiperiodic case as well, and indeed the 
numerical results show the expected scaling. This al- 
lows one to extract the infinite size limit of the ground 



QMC 
gs,OT 



-0.6581(1), to be compared 



onal tiling [46[ is E'^ 

with the square lattice value (4^ of -Eg^ff = -0.66943. 
For the Penrose tiling, the GSE found by QMC [3i is 



E 



QMC 
gs,PT 



-0.6529(1). All energies are expressed in units 
of J, the exchange coupling. One sees that the ground 
state energy of the square lattice lies below that of the oc- 
tagonal tiling which lies below that of the Penrose tiling. 
As the classical values for all three systems is exactly 
equal, Eg\ = —0.5, the differences are due to the dif- 
ference in quantum fluctuations in these three types of 
antiferromagnets. This point will be discussed again in 
the section devoted to the energy spectrum in linear spin 
wave theory. 

The staggered magnetizations cannot be directly mea- 
sured, as one is dealing with finite systems, for which 
there is no symmetry breaking. The expectation value 
{Siz) giving the staggered moment on site i is zero. In 
periodic lattices, it is nevertheless possible to obtain esti- 
mates of the values of the local staggered moment at each 
site by summing the pair correlation function over all 
pairs of spins. This estimator becomes exact in the ther- 
modynamic limit, N — > 00. Generalizing to the tiling, 
one defines 



QMC 



3 ^ 



(5) 



where the sum extends over all sites and the SU{2) 
symmetry of the Hamiltonian was used. The quantity 
^QMC (^j^ is the QMC estimate of the staggered magne- 
tization at lattice sites i. Measurements of this quantity 
for the octagonal tiling are plotted against the coordina- 
tion number z in FiglT] (data obtained for the N=1393 
approximant). The values shown were normalized with 
respect to the average value Jl2 — 0.337 ± 0.002 deter- 
mined by extrapolation to the thermodynamic limit. 

As seen in FiglJl local environments clearly play an 
important role. To sort out the roles of short-range versus 
long-range effects, one can introduce another estimator 
for the local staggered magnetization. 



QMC,-. _ 




(6) 



where the sum is truncated to the sites in the immediate 
vicinity of i. This quantity is a local approximation of 
the one defined in EqlS] The results for this local estima- 
tor for the local magnetization is plotted in Fig|S] along 
with, for comparison, the prediction of the Heisenberg 
star calculation (Eq|4]). The trend towards smaller values 
of staggered magnetization as z is seen in both cases. One 
can conclude from these results that the nonmonotonic 
dependence of the full staggered magnetization M;^c^'^(*) 
as a function of z is due to longer range correlations in the 
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FIG. 7; The normalized staggered moments plotted against 
coordination number z (calculation for N=1393 approximant 
of the OT). The inset shows the spatial variation of the stag- 
gered magnetization where the radii of the circles correspond 
to the size of the local staggered magnetizations. (Reproduced 
from gi]) 



quasicrystal. This will be further explored in the section 
devoted to linear spin wave theory. 




0.78 



0.73 



0.68 



0.63 



0.58 



FIG. 8: The local magnetization plotted against coordination 
number z. The points on the dashed line are the exact values 
obtained for Heisenberg stars. The inset shows a portion of 
the tiling where the radii of the circles correspond to the size 
of the local magnetizations. (Reproduced from [4^) 

The global average of the staggered magnetizations 
for the octagonal tiling is /is = 0.337 ± 0.002 which 
is larger than the value found for the square lattice 
{fis = 0.3071 ± 0.0003[44]). This indicates that quan- 
tum fluctuations reducing the order parameter are sup- 
pressed in the quasiperiodic tiling due to the inhomoge- 
neous connectivity. The ground state energy per spin is 
—0.6581, compared to the value —0.6694 (square lattice) 
again showing that the square lattice has bigger quantum 



fluctuations than the quasicrystal. 

For the Penrose tiling, QMC calculations [i^ for the 
staggered moments show that, as for the octagonal tiling, 
staggered magnetizations are distributed in a range of 
values, with a strong coordination number dependence. 
These results are reviewed later on, when we compare the 
QMC and linear spin wave theory for the Penrose tiling 
antiferromagnet . 

To explain the observed results one can turn to other 
theoretical tools such as renormalization group or linear 
spin wave theory, described next. 



C. Scale transformations and Renormalization 
Group 

We have mentioned that an important symmetry of the 
tiling, namely invariance under discrete scale transfor- 
mations (inflations/deflations) suggests the possibility of 
calculating properties by a renormalization group (RG). 
This method is used for problems with scale invariance, 
and has been successfully applied to systems near a crit- 
ical point, for example. 

The renormalization group scheme described here is a 
generalization of the calculation of Sierra and Martin- 
Delgado for the square lattice antiferromagnet [i^. In 
this approach 32, l33| , it was shown that under a change 
of scale, the original Hamiltonian is transformed into a 
new Hamiltonian of the same form but with rescaled pa- 
rameters. The first step consists of defining blocks of 5 
spins (formed by a central spin S and its four nearest 
neighbors), such that the blocks form a square lattice 
of length that is ^/E times bigger. The original spins 
of each block were next combined to form an effective 
block spin. The important point is that the interactions 
between the block spins can be written in a form of a 
Heisenberg model with a new nearest neighbor coupling 
which depends on the old coupling, as well as the spin and 
length rescaling factors. The mapping is not exact, and 
among other approximations, some terms of the Hamilto- 
nian are neglected. This calculation for the square lattice 
was adapted to the octagonal tiling in [sl]. The case of 
the Penrose tiling, more complex, has not so far been 
attempted. For the octagonal tiling, the global ground 
state energy was computed, and in addition, values for lo- 
cal site-dependent order parameters were obtained. The 
approach is briefly described below. 

We have seen that the sites of the octagonal tiling are 
classified into A,B,C,D,E and F in order of decreasing 
coordination number (from z = 8 for A sites down to 
z = 3 for F sites). The frequencies (or probability of 
occurrence in the tiling) of each type of site is exactly 
known, in terms of the number A. Consider a finite piece 
of the tiling of N sites. Under deflation, where the length 
scale is increased by A, the number of sites is reduced by a 
factor 1/ A^ . One can show that the sites which disappear 
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are those having z = 3,A and half of the sites with z = 5, 
called D2 sites. The sites that remain are the A,B,C, and 
half of the D sites (the Di sites). They are transformed 
into the A,B,C,...,F sites of the new tiling and have a new 
value of the coordination number z' < z. For example, a 
B site always transforms into a I?2 site under deflation, a 
C site site into an E site, and so on as in Table 2. These 
four families of sites are used to introduce block spins 
variables. 



initial site 


final site 


(z value ) 


(z' value) 


A (8) 


A,B,C or Di (8,7,6,5) 


B(7) 


-D2 (5) 


C(6) 


E (4) 


^1 (5) 


^ F(3) 



Table 1. Transformations of sites in the octagonal tiling 
under an inflation 



Each of the "new" sites after deflation is associated 
with a given cluster of sites, comprising a central site and 
its z neighbors, on the original tiling. There is however 
a first complication that was not present in the square 
lattice, due to the fact that, in the octagonal tiling, some 
of the clusters turn out to share edges. The first approx- 
imation made in the RG scheme is therefore to remove 
a small fraction, ~ 10% of the bonds. This bond 

dilution allows decoupling of clusters. Fig|9] shows the 
sites of the block spins by black dots, and the bonds 
that have to be taken out in order to obtain decoupled 
clusters. Note that the dilution of the OT leaves it two- 
dimensionality intact. 




FIG. 9: Tiling showing block centers (black dots). The grey 
lines connect pairs of sites that are shared between two blocks. 
(Reproduced from [sl]) 

One can now define a new "block spin" variable S' = 
{z—V)S composed of the old spins S of the z -I- 1 spins be- 
longing in each of the clusters. FigfTOl shows an F-site (of 
coordination number z = 3) on the original tiling, and, 
on a larger length scale, a composite F-site made up of a 
central block spin surrounded by three other block spins 
on the tiling. Upon repeating the infiation transforma- 
tion, clusters can be defined on larger and larger length 



scales. Notice that each of the four block spins of the 
figure on the right has a different internal structure - 
this illustrates one of the differences that make the qua- 
sicrystal such a difficult problem compared to a periodic 
system. The next important approximation made in the 
RG is to restrict the number of such combinations possi- 
ble, so that one works with a finite set of variables rather 
than the infinite number that is needed in principle for 
the quasicrystal. In the Heisenberg model, this trunca- 
tion to a finite set of environments works because of the 
fact that next-neighbor and further-neighbor corrections 
decrease very quickly on the tiling. Contributions to the 
effective Hamiltonian of a given central spin due to far-off 
spins are quantitatively unimportant, and they will sim- 
ply lead to a broadening of the block spin values obtained 
for each of the five types of blocks. 




FIG. 10: (left) Example of an F-site(2=3)(left) with the cen- 
tral and peripheral sites indicated by red dots, (right) Exam- 
ple of an F-site on a larger scale with block spins defined on 
clusters of sites as shown. (Reproduced from [33| ) 

As in the square lattice, it can be shown that the block 
spins of the inflated tiling interact via an effective an- 
tiferromagnetic interaction. Just as in the case of the 
square lattice [i^, one can write an approximate effec- 
tive Hamiltonian as a sum of single-block contributions 
(the diagonal terms) and a set of terms involving near- 
est neighbor blocks (off-diagonal terms). That is to say, 
the expression for the energy of the inflated system can 
be written so as resemble the original nearest neighbor 
Heisenberg Hamiltonian, up to multiplicative factors and 
an additional term. The renormalization group (RG) 
transformation relates the new and the old Heisenberg 
Hamiltonians, by specifying a) the set of transformations 
of spins S into block spins S'^ (where j labels the type 
of block) and b) the set of transformations of the old 
couplings {J} into new effective couplings {J^j}. The 
number of spins is reduced from N to TV'. The formal 
expression for the transformed problem reads 

tIh(N, S, J)To « iV'eo(J, S) + H'{N', S', J') 

(7) 

where, after the scale transformation Tq, the new 
Hamiltonian H' has the same form (bilinear in S') as 
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H. The new problem that one encounters in a qua- 
sicrystal, as we have seen already, is that the number 
of types of blocks and the number of types of couplings 
between them increases with the number of steps of in- 
flation. The approximate scheme described in [33| trun- 
cates the number of cluster types and cluster-couplings 
that are retained in the calculation to just seven cluster 
types and five inter-cluster couplings. The new couplings 
couplings are related by a multiplicative (matrix) re- 
lation to the old and and decays with increasing number 
of deflations. The block spins are similarly related by a 
matrix equation to the original spins and they increase 
with n. The resulting set of matrix equations for the 
couplings and the spins is then diagonalized to obtain 
the fixed point solution of the system. The solutions for 
the block spins flow to infinity, i.e. the classical limit, 
as in the periodic lattice, with a fixed ratio between the 
seven different block spins. 

The ground state energy of a single block is that of 
a z-fold Heisenberg star (Sec. II.??) in which z spins of 
lengths Sj are coupled with strength J to a central spin 
Sq. The total ground state energy Eq can be written as 
the sum over all blocks of the block energies, at all orders, 
as follows: 



+ ) (8) 



The zeroth order contribution is the sum over cluster 
energies ef^ where i labels the type of site on the original 
tiling, whose frequency of occurrence is N'^^^{z)/N = fi. 
With each deflation, the number of clusters decreases as 
A^. The energies after n deflations, e^"'\ are the energies 

(n) 

of blocks with a spin Sq at the center, with effective 
couplings J^"^ to the S'|"'' surrounding spins. Summing 
the series for the energy gives Eq/N sa —0.51, which lies 
about 20% below the value deduced from numerical QMC 
data fis'l . This order of magnitude of error is also found 
in the RG solution for the square lattice, and arises in 
part due to the neglect of certain types of bonds in this 
approximate scheme (see (igj). 

For the tiling, it is possible to calculate the relative 
strength of the local staggered magnetizations on each of 
the six families of sites by this RG method. These are 
determined by converting the cluster energies into esti- 
mates for the staggered magnetization using the same 
argument as in EqH) Fig (TT] shows the staggered magne- 
tizations m.i (as a fraction of the full value S) for z values 
ranging from 3 to 7 for a 4th order calculation after cor- 
recting for the bond dilution. Also shown are the zero 
order result (dashed line) and the QMC data. The order 
parameter decreases with increasing z, agreeing qualita- 
tively with Quantum Monte Carlo . The quantitative 
accord is also very good. 
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FIG. 11: Staggered magnetization (as a fraction of the full 
value 5") obtained by the RG calculation rrii at fourth order 
after correction for bond dilution (grey rectangles), the QMC 
data (grey circles)and the analytic zeroth order result (dashed 
line). (Reproduced from [3^ ) 



D. The Linear Spin- Wave Approximation 

Spin wave theory gives the quantum corrections to 
classical ground state. It is only useful in situations when 
the exact ground state is known to be close, in a sense 
to be clarified below, to the classical state, in this case, 
the Neel two-sublattice up-down state. In the quasicrys- 
tal, as in other unfrustrated Heisenberg antiferromag- 
nets, one expects there will be long range order at zero 
temperature, as in the previously described sections. In 
the ground state, spins acquire a non-zero average value 
along the z-direction, with mosi — {eiSf)o ^ where 
Ci = ±1 depending on whether the site is on the A or 
B sublattice. In spin wave theory, one assumes that the 
deviations from the classical value are small compared to 
the spin length, i.e. S — mosi << S, justifying a pertur- 
bation theory around the classical state. In this case the 
exact Heisenberg Hamiltonian can be approximated by 
the first several terms of a power law expansion in spin 
deviation operators. Quantities such as magnetizations 
and bond strengths can then be calculated to successive 
orders of this expansion. 

In the quasicrystal, the sites are expected to have dif- 
ferent values of the local staggered magnetization. We 
assume that quantum fluctuations are small for each of 
the sites taken separately, and carry out an expansion 
around the classical antiferromagnetic ground state by 
introducing spin deviation operators for each site. Ac- 
cording to standard definitions (see for example (40|). 
this means introducing two families of Holstein-Primakoff 
boson operators: ai,a| for sites on the A sublattice, 
bi,bl for sites on the B sublattice. The boson operators 
account for deviations from the classical configuration, 
with the spin deviation operator on A-sublattice sites 



given by a|a 



rii = S — SI and on B sublattice sites 
S + Sj . Expanding the original Heisenberg 
Hamiltonian and keeping terms upto second order in the 
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boson operators one obtains the linear spin wave (LSW) 
Hamiltonian 

Hlsw = -JS{S+l)Nb+JS (ala^+bjbl+albl+bja^) 

<ij> 

where N\j — 2N is the number of bonds, with N being 
the number of sites. This quadratic Hamiltonian can be 
diagonalized by a generalized Bogoliubov transformation 
[50| . The diagonalization is carried out in the real space 
basis, for finite systems ~ contrarily to periodic systems, 
where using a Fourier representation leads to significant 
simplification of calculations. 

After the Bogoliubov transformation, which takes the 
set {tti^bj } to a set {am,Pm }, the Hamiltonian takes 
the diagonal form : 

N/2 

Hlsw = ^0+^ n-^al^a^ + n:^^Pl^(3m (9) 

m— 1 

where and fl^ are the energies of the eigenmodes. In 
crystals the eigenmodes described by the /3 operators are 
spin waves: they carry both a spin and a crystal momen- 
tum. In the quasicrystal, the eigenmodes are spatially 
much more complex. Depending on the energy of the 
mode, we will see that they can be more or less spatially 
extended, and for certain values of the energy, they can 
be strongly localized on a subset of sites of the tiling. 
Nevertheless, in the limit of very low energies, or very 
large length scales, the details of the quasiperiodic struc- 
ture are not pertinent, and the solutions are expected to 
approach extended states. In this limit, the dispersion 
relation should resemble the usual linear relation found 
in antiferromagnets when /c — >■ 0, E{k) « cfc, along with 
a density of states N{E) ^ E'^ in d — 2. This is indeed 
observed (see the section on the density of states) . 

Eq is the ground state energy, comprising a classical 
part and a term arising from quantum corrections SE, 

N/2 

EI^^^ = -JS{S+l)Nb + JSY,^nr 

m— 1 

= -JS^Nb + SE (10) 

The local staggered magnetizations are calculated using 
the following expressions valid for sites i{j) belonging to 
the A(B)-sublattice respectively: 

ms{z)^\{S^)\=S-{ala,)o 

m^ij) ^ \{S;)\ = S - {b]b,)o (11) 

where the expectation values are calculated in the ground 
state. 

E. Finite size samples for numerical calculations. 

Calculations can be done for finite samples of two 
types. One can consider finite portions of an infinite sys- 



tem, with free boundary conditions. Such samples have 
the usual problem of spurious boundary modes, which 
are difficult to correct for. Such samples can be gen- 
erated by projecting down from five dimensions (in the 
Penrose tiling) or four dimensions (octagonal tiling) us- 
ing the standard projection matrices. The second alter- 
native is to consider approximants as in previous work 
on electronic (see review [sH) or phonon models [sl ]. 
These, as we saw already are closely related to the in- 
finite quasiperiodic structure at small length scales, but 
different on large length scales. For the Penrose tiling, 
one can generate rectangular shaped approximants obey- 
ing periodic boundary conditions called the Taylor ap- 
proximants For the octagonal tiling, one can gen- 
erate square periodic approximants [s^. By considering 
the problem in periodic approximants one eliminates the 
problem of boundary effects, at the cost of introducing a 
global shear-type distortion. The latter has consequences 
for the energy spectrum as the approximant has a sym- 
metry not present in the quasicrystal. It turns out that, 
fortunately, there is good agreement between the results 
for the two different kinds of system when one consid- 
ers thermodynamic quantities obtained by averaging over 
modes. For example, the order parameters for spins deep 
in the interior of finite samples have the same distribution 
as in the periodic approximants. 



V. RESULTS FOR GROUND STATE ENERGY 
AND STAGGERED MAGNETIZATIONS 

For the Penrose tiling, the ground state energy per site 
(GSE/site) is found from EqlTU]for the periodic approx- 
imants of size N ranging from 96 to 11556. The extrap- 
olation to infinite size gives an asymptotic value Eq — 
-0.643(0) ± 0.0001. For the octagonal tiling, the LSW 
calculation gives a ground state energy Eq — —0.646 |53| . 
The value for the square lattice is -0.658 (which decreases 
to -0.6705 when the next term in the perturbation series 
is taken into account, see (40|'). 

The ground state energy per site and bond is shown 
in Table 1 for different types of lattices. Of the three 
systems presented here, the lowest GSE/site is that of 
the square lattice, then the octagonal tiling, and, very 
slightly higher, the Penrose tiling. In addition, we show 
an approximant of the octagonal tiling (the "crown" lat- 
tice), which has 7 sites per unit cell and was studied in 
[H^ . This approximant intercalates between the square 
lattcie and the octagonal tiling. The ground state of 
the Penrose tiling is closest to the classical ground state 
because, presumably, spin waves are "harder to excite". 
These results should be treated with caution because 
they are obtained for approximant systems. The higher- 
order corrections in spin wave theory will play an impor- 
tant role in dtermining the final value of ther spin wave 
velocities. Nonetheless, based on the LSW results, the 
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higher energy of the PT could reflect the slightly higher 
"complexity" of the PT compared to the OT. One could 
expect, therefore, that the GSE/bond for tilings of high 
co-dimensions [38[ will lie above that of the PT. The 
GSE/bond of the Rauzy tiling of Vidal and Mosseri [sl], 
would, on the contrary, lie below that of the OT. These 
expectations have not been checked so far by calculations. 



lattice 


GSE 


GSE 


bond 


GSE/bond 




(cl.) 


(LSW) 


(cl.) 


(QMC) 


Penrose 


-0.5 


-0.643 


-0.25 


-0.322 


Octagonal tiling 


-0.5 


-0.646 


-0.25 


-0.323 


Crown lattice 


-0.5 


-0.648 


-0.324 


-0.330 


Square lattice 


-0.5 


-0.658 


-0.25 


-0.335 



Table 1. Classical and LSW values for ground state energy 
per site and per bond (for S = 1/2 and in units of J) for the 
two quasiperiodic antiferromagnets, an approximant system 
(the crown lattice) and the square lattice 

Fig HH indicates how values of the local magnetizations 
vary with z value, in the Penrose system (Taylor approx- 
imant size N — 4414). The figure shows the average val- 
ues and the standard deviations of mi{z) for each value 
of 2, as obtained in LSW theory and from QMC , as re- 
ported in [H^. The figure shows that Quantum Monte 
Carlo gives a narrower spread of the magnetization. Lin- 
ear spin wave theory clearly overestimates the fluctua- 
tions from average behavior, giving a too high value for 
z = 3 and too low values for higher values of z. The 
LSW theory tends to exaggerate the spread in values of 
the staggered magnetizations compared to QMC. FigfT^ 
shows the details of the rrii (z) distribution for each of the 
z values. There is a continuum of values, as expected for 
a quasiperiodic structure, with some pronounced peaks. 
Each of the main peaks correspond to specific next near- 
est configurations, as we will discuss next. 

In order to explain the multiple peaks in the values of 
nisi seen numerically in Fig |13l for a given z, we must take 
into account medium and long range structural details. 
This can be done analytically by considering clusters like 
the Heisenberg star cluster introduced already. As was 
shown in [H^] , it is easy to carry out a LSW type analysis 
for the simple z-leg Heisenberg star in FigEl and calcu- 
late the magnetization of the central site as a function 
of z. To understand the results for the Penrose tiling, 
one needs to go a step further, and consider the effect of 
next nearest neighbors. This can be done by considering 
a two-tier Heisenberg star, as shown in FigHH The LSW 
result obtained in 56i] showed that the staggered magne- 
tization which now depends on z and z', ms{z, z') has a 
nonmonotonic behavior in z. The fact that the distribu- 
tion of TTis for the z — 3 sites displays three peaks can 
be explained being due to the three different categories 
of these sites when one takes next nearest neighbors into 
account. Similarly, the z = 5 sites fall into two categories 
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FIG. 12: Averages and standard deviations of the local 
staggered magnetization as a function of z ( LSW theory 
(dashed line) QMC (continuous line) for Taylor approximant 
N=4414). (Reproduced from [H) 
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FIG. 13: Probability distribution of the local staggered mag- 
netizations for the five different values of z in a Taylor ap- 
proximant (N=11556). (Reproduced from fsst]) 



due to the difference in z', resulting in the splitting of the 
distributions into two groups for this class of sites. 

Figim shows by an intensity plot how local magneti- 
zations vary in space on a portion of the Penrose tiling. 
The lowest values of staggered magnetization are found 
on a certain subset of z = 5 sites, while the largest values 
are found on the sites of z = 3. 

A. Self symmetry of the staggered magnetization 

We have seen that the values of rusi for a given z are 
distributed about a mean value that depends on z, and 
that the actual shape of the distribution depends on de- 
tails of the next nearest neighbors, and further nearest 
neighbors. In this section, we focus on the order pa- 
rameters of a subset of sites which are self-similar under 
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FIG. 14: Two-tier Heisenberg star showing z near neighbors 
and z' next nearest neighbors of the central site. (Reproduced 
from [H) 



FIG. 16: Self-similar structure of the order parameter: the 
y-axis shows niai values of 40 A-sites of the octagonal tiling 
numbered (x-axis) according to increasing niai value. The 
values of the coordination number after one inflation {z') and 
two inflations {z") have been indicated. (Reproduced from 




FIG. 15: Intensity plot of the staggered local magnetization 
nis in the Penrose antiferromagnet. Lighter and darker shades 
correspond to local values which are greater or smaller than 
the average value, respectively. 



inflation and we show the resulting self-similar structure 
of the distribution of order parameters. These are the A- 
sites of the octagonal tiling, i.e. the sites with eight-fold 
symmetry and coordination number z = 8. Calculations 
show that the nisi values for A-sites are grouped into four 
disjoint groups, each of which corresponds to a differ- 
ent next nearest environment, as we already explained. 
The uppermost group of values corresponds to A-sites 
which transform under deflation into A-sites. The next 
lowest group of nisi values correspond to A-sites which 
transform into B-sites (z' = 7), etc. Now, as can be 
seen in FigHni the uppermost group can be seen itself to 
be composed of four subgroups, of which the uppermost 
subgroup corresponds to A-sites which transform under 
deflation into A-sites. This substructure would continues 
infinitely in an infinite system but stops at the third level 
of the hierarchy in the small samples that were considered 
for the numerics. 



B. Perpendicular space representation 

It is interesting to consider the staggered magnetiza- 
tion from the viewpoint of the perpendicular space of 
the tiling [s^. Recall that, in perpendicular space, the 
projections of different families of sites fall into separate 
domains. When the msi values are plotted as a function 
of fj_ rather than r^j ( as was done in Fig llSI) one obtains 
the figures shown in Fig|17l The four domains shown 
lie at different distances along a third axis (not shown). 
Instead of the complex variations see in the real space 
representation, one sees five-fold star shaped domains, as 
illustrated in FigllTl The figures show that, as expected, 
the domain corresponding to each coordination number 
has a distinct color. The colors are not absolutely uni- 
form since no two sites are identical. 

Self-similar patterns can be investigated using this type 
of representation. To illustrate how this occurs, consider 
the octagonal tiling, whose behavior under infiation was 
discussed in the section on RG. If one considers the A 
(eight-fold) sites, for example, their transformation under 
inflation depends on the position in perpendicular space. 
The closer r± to the origin, the bigger the number of in- 
flations under which A ^ A. Another way of saying this 
is that for such sites (i.e. those close to the origin of the 
perpendicular space window) the environment remains 
eight-fold symmetric out to greater distances. The stag- 
gered magnetization which depends on the environment, 
thus, depends on the perpendicular space coordinate, and 
should vary in magnitude progressively as one goes out- 
ward from the origin. To measure accurately this self 
similarity, flrst described in [46], it requires calculations 
on systems large enough to permit a large number of in- 
flations and to overcome the effects of distortions due to 
the boundary conditions. 

To resume, the intensity plots in perpendicular space 
can serve not only to represent in an accessible and com- 
pact form the complex real space distribution of the 
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order parameters on the quasiperiodic tiling, but also 
to demonstrate the properties of the state under length 
rescaling. 



3f 



FIG. 17: The same data as in Fig llSI for the ground state 
staggered magnetizations (ms) in perpendicular space. Each 
of the vertices of the Penrose tihng project into one of four 
planes (see Box 3). The intensity at the point represents the 
value of ms on that vertex. 



Static structure factor 



The static longitudinal structure factor is defined by 



N 



-ik.(fi —fj ) 



(SziSzj) 



(12) 



Figim shows an intensity plot of this static structure 
factor calculated for an approximant of the Penrose struc- 
ture, for the nonmagnetic state and the antiferromag- 
netic ground state. There are many peaks of very low 
intensity of S. The figure shows only the most intense 
peaks, whose positions are indicated by circles (squares) 
for the magnetic (nonmagnetic) structures respectively. 
The five-fold symmetry present in the magnetic state, as 
it is in the nonmagnetic state, along with the symmetry 
under fc — )■ — fc lead to the observed 10-fold symmetry of 
the peaks. The octagonal tiling antiferromagnet has a 
structure factor with 8-fold symmetry of both the mag- 
netic and the nonmagnetic states. In both tilings, it can 
be seen that the magnetic peaks lie between the nonmag- 
netic peaks with no coincidences between the two. Peaks 
can be indexed using sets of four (or five) indices. These 
indices have integer values for non-magnetic peaks and 
can be integer or half-integer for the magnetic peaks. As 
in the case of the octagonal tiling [s^ the half integer 
indices arise due to the doubling in size of the antiferro- 
magnetic unit cell in five dimensions. This, coupled with 
the extinctions in the structure factor rule due to the mul- 
tiplicity of each unit cell, then leads |53| to the observed 
"shifting" or "displacement" of the magnetic peaks with 
respect to the nonmagnetic ones. This shifting of indices 
would be true also for a three dimensional bipartite sys- 
tem, such as the icosahedral tiling, with spins interacting 
along edges of the rhombohedra. If all of the sites were to 
be occupied by spins, the antiferromagnetic state would 
give rise to half-integer indices as well in the structure 
factor. However, real life compounds such as ZnMgHo 
are more complicated. This compound was studied by 



FIG. 18: Intensity plot of the static longitudinal magnetic 
structure factor ^"(k) for the — 4414 Taylor approximant. 
The highest intensity peaks are shown in the figure at posi- 
tions that are indicated by circles (squares) for the magnetic 
(nonmagnetic) structures respectively. The relative intensity 
is denoted by a linear gray scale ranging between zero (white) 
and maximum intensity (black). (Reproduced from [551 ]') 



neutron scattering [14], and seen to have short range 
antiferromagnetic correlations below about 20K. These 
correlations lead to a magnetic superstructure that is, 
as for our two dimensional model, shifted with respect 
to the paramagnetic state. The antiferromagnetic vector 
that best fits the data has a more complicated value than 
the simplest form for a 3d quasiperiodic antiferromagnet 
{q,^^,i = 1,6). This may occur because the magnetic 
unit cell is much larger for the three-component system, 
due to the fact that only the Ho sites carry a magnetic 
moment, and this results in smaller - fractional - dis- 
tances between peaks in reciprocal space. 



VI. MAGNON SPECTRUM AND 
WAVEFUNCTIONS 

A. Energy spectrum and density of states. Linear 
dispersion law. 

The low energy spin excitations, or magnons, have en- 
ergies f2^. The two sets of energies f2+ and fi" become 
identical in the limit of infinite size when sublattices A 
and B become strictly equivalent. The integrated density 
of states (IDOS) for the two sets of modes is defined by 



1 ^ 



(13) 



m— 1 



where a = ±. The IDOS is plotted in FigHlfor three 
different system sizes. The salient features of the figure 
are the several groups of closely spaced energy levels, the 
main gaps, which are stable with increasing system size, 
and a plateau at the energy E = 3, corresponding to a 
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FIG. 19: a) (Color online) Integrated density of states N{E) 
versus E (expressed in units of J) calculated for three con- 
secutive approximants ( black dots : A*" = 644, red dots : 
N = 1686, grey : = 4414) . b) Low energy tail of N{E) 
versus (the straight line is a fit to the data). (Reproduced 
from HI) 




FIG. 20: Density of states calculated for the Taylor 6 approx- 
imant (N=11556 sites). (Reproduced from [sF 



the 5-fold sites. 3) The peak at the exact value oi E = 3 
corresponds to string-like states living on the z = 3 sites. 
4) The lowest energies correspond to the most extended 
states of the spectrum, having a two dimensional charac- 
ter and having its support on all types of sites. 

As shown in Figll9b) the low energy part of the IDOS 
can be fitted to a power law of the energy, N{E) = 
(sttc^) This is consistent with a linear dispersion of 
the magnon modes in this region of the spectrum. Fitting 
to a form N{E) — E"^ /{%ttc^) gives a spin wave velocity 
on the Penrose tiling of c = I.IJ. It is interesting to 
compare this result for c with the corresponding values 
on the square lattice and the octagonal tiling. On the oc- 
tagonal tiling, the value is Cocta ~ 1-3 J (Nb. a lower value 
of 1.1 J was obtained from a fit to the dynamical struc- 
ture factor 53|). The error in the spin wave velocities of 
the quasiperiodic systems is large, of the order of 10%, as 
the value depends on the region of the fit, and the system 
sizes that were considered are not big enough to allow a 
good smoothing of the low energy part of the IDOS. On 
the square lattice, the value can be calculated analyti- 
cally in LSW theory, and it is Cgq = 2V2JSa « 1.41J 
(for edge length a = 1 and S — 1/2). 

To resume, excitations in the quasiperiodic antifcrro- 
magnets obey an acoustic-type dispersion relation at long 
wavelengths, as expected on general grounds for antifer- 
romagnets. In LSW theory, the values of the spin wave 
velocity in the Penrose tiling is smaller than in the oc- 
tagonal tiling, which is in turn smaller than the value 
on the square lattice. For all systems, the edge lengths 
are all taken to be equal to unity. The propagation of 
modes appears to be slowest in the Penrose tiling. For 
the same edge length, one can calculate that the density 
of vertices p (number of sites per unit area) is greatest for 
the Penrose tiling (p=1.2311), followed by that of the oc- 
tagonal tiling (/3= 1.2071), followed by the square lattice 
(/9=1). Since at long wavelengths the spin velocity de- 
pends on the mean density, this provides an explanation 
of the trend observed for the sound velocities in these 
three systems. 



B. Wavefunctions. 



set of exactly degenerate states which will be discussed 
in the section on wavefunctions. 

The derivative of the IDOS is shown in Figl20l One 
can distinguish several groups of energies separated by 
gaps. The local density of states (LDOS) allows one to 
determine which sites are most important at any given 
energy E. A study of the LDOS for the Penrose tiling 
system was carried out in [5^ and leads to the following 
conclusions : l)The two highest energy bands around the 
values E « 5.4 and E « 6.4, correspond to wavefunctions 
that are localized on the 6- and 7-fold sites. 2) Lying be- 
low these in energy are the states centered primarily on 



We now discuss the magnon wavefunctions and their 
spatial characteristics for the different parts of the energy 
spectrum. 

As seen by the large step in dN{E)/dE at E = 3 
(Fig|20|). a large number of degenerate states occur at 
this energy. This is due to wavefunctions that have their 
entire support on the z ~ 3 sites. These are string-like 
states forming closed loops. The smallest loops are rings 
around the "football" z = 5 sites. One typical linear com- 
bination of such degenerate states is shown in Figi2Tb). 

The band of highest energies corresponds to states that 
have large amplitudes on sites of z = 7. The highest en- 
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FIG. 22; The inverse participation ratio for states in OT ap- 
proximants plotted as a function of energy (uj/ JS) for three 
different system sizes. Inset: the IPR as a function of en- 
ergy for square lattice samples of different sizes, as indicated. 
(Figure reproduced from [53l ]) 




FIG. 21: Intensity plots in a PT approximant (N=4414) show- 
ing the probability ipi{E)^ for states of energy E (a darker 
shade corresponds to a higher probability): a) _E = 6.469 b) 
E = 3.0000 c) 0.8887 d) E = 0.1692. (Reproduced from 
M) 



ergy band has a width of about 0.16 and is centred around 
E = 6.4. The states involving z — 6 sites correspond to 
energies in the range 5.43 < E < 5.49. These states 
have a smaller dispersion (flatter band) than the states 
in the topmost group of energies. This may be explained 
by noting that the six-fold sites are significantly fewer in 
number than the seven-fold sites. These sites are there- 
fore spaced further apart, which should lead to smaller 
wavefunction overlaps and smaller bandwidth. 

The lowest energy states , for i? < 1, correspond to 
wavefunctions for which the amplitude is non-zero on 
all categories of sites. The wavefunctions are in this 
sense more two dimensional than the previously men- 
tioned wavefunctions. On the other hand, they are not 
truly extended - they have large amplitudes within re- 
gions of size that depend on the energy i?, and fall rapidly 
to zero outside these. It would be interesting to carry de- 
tailed studies of the behavior of these low energy states 
for large systems, in order to explore their multifractal 
properties. Figsi2Th) and d) show the states correspond- 
ing to two energies. 

In disordered systems, a quantity that is often com- 
puted in order to determine the localized or extended 
character of wavefunctions is the inverse participation ra- 
tio (IPR), defined by 



P-\E) 



The inverse of this quantity, or PR, corresponds, 
roughly speaking, to the number of sites on which the 
wavefunction is finite. For an extended state, this num- 
ber is of the order of the system size N, while for a lo- 
calized state the PR is a much smaller number - only 
a few sites participate in the wave function. In general, 
the scaling of the IPR with system size is a gauge of 
the degree of localization (this is not true when there 
are macroscopic degeneracies, as in the string states of 
E = 3). If the IPR varies as P^^ « N'l^, then, as 
seen from the preceding remarks, (3 ~ \ corresponds to 
extended states, /3 = to localized states and intermedi- 
ate values of beta to other cases including the so-called 
critical states. 

Milat and Wessel have calculated the IPR on three ap- 
proximants of the octagonal tiling, and these are plotted 
as a function of the energy in Figl22l One sees that, at 
low energies, the IPR depends on the system size, an in- 
dication that the corresponding states are extended. At 
high energies, the IPR flattens out and less size depen- 
dent, indicating that states are more localized. The inset 
shows, for comparison, the IPR calculated for square lat- 
tice samples, in which one knows the eigenstates to be 
Bloch states and therefore perfectly extended. 

The studies of the IPR on approximants of the Penrose 
tiling 155)] confirm the observation that states are more 
delocalized at low energy compared to higher energies. 
The value of /3 are largest in the low energy end of the 
spectrum, where it has an average value of slightly less 
than 1 (/3 = 0.9±0.1). The multifractal properties of the 
eigenstates can be investigated by calculating the func- 
tion /(a) which contains information on the scaling of 
the wavefunction - f{a) shrinks to a single point for a 
simple fractal or completely extended state but has a typ- 
ical convex shape for Cantor set like wavefunctions, such 
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as those at the metal- insulator transition [58|- There is 
some evidence for a nontrivial /(a) spectrum [55|, and 



thus multifractal states, in the low energy region of the 
spectrum. 

The figures showing wavefunctions in real space can all 
be converted to perpendicular space representations [H^ , 
where there is a change, from localized to delocalized, of 
states as a function of E. This provides an alternative 
way of probing the wavefunction properties. 

To sum up, we have shown some representative exam- 
ples of magnon wavefunctions for the Penrose tiling ap- 
proximants. Magnon wavefunctions have not been stud- 
ied in sufficient detail, in large enough samples, how- 
ever. Another problem that we have already mentioned 
in the context of the energy spectra is that relating to 
the boundary conditions. The wavefunction calculations 
shown here were done on a Taylor approximant of 4414 
sites, with periodic boundary conditions. It is clear that 
the boundary conditions will have significant effects on 
the wavefunctions, particularly at small energies, where 
the wavefunctions are the most extended, although all 
other eigenmodes are in principle affected as well. This 
important problem has not been addressed satisfactorily 
thus far in any of the numerical calculations on tilings, 
be it for spins or electronic models. 



C. Dynamical structure factor 



The dynamical magnetic structure factor is defined by 
1 ^ 



di(S^,(t).S^,(0)) (14) 



where S_Li = Sxi, Syi denotes the transverse components 
of the spin Si. For w = 0, the dynamical structure fac- 
tor reproduces the magnetic peak structure of the static 
structure factor already defined in Eq[T21 The u!{k) can 
be measured by scans in fc-space. FigsHS] show the in- 
tensity of the structure factor for scans effected around 
different magnetic Bragg peaks. The figure on the right 
shows the intensity profile in the k, uj plane as one varies k 
along the line joining two intense magnetic Bragg peaks 
at 1.27r, 0.57r and 37r,37r. One can see the soft modes, 
and linear dispersion near these two peaks. A fit to the 
data gives a spin wave velocity c of about 2.1, in accord 
with the value that was extracted from the fit to the low 
energy end of the IDOS. 



D. Local frequency-dependent susceptibilities 

The dynamical structure factor yields information on 
the imaginary part of the local susceptibility via the re- 
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FIG. 23: Intensity plot of the dynamical structure factor of an 
OT approximant (N=1393) in the vicinity of magnetic Bragg 
peaks (left) along the line joining the peaks at fcs — 1.2n, O.Stt 
and ks' = 37r, 37r and (right) along a line parallel to the ky 
axis and traversing ks = 2.9it, 4.87r(right) (Figure reproduced 
from HI) 
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FIG. 24: Imaginary part of the local dynamic susceptibility 
plotted as a function of frequency lu averaged separately over 
sites belonging to six different z values on the OT (Figure 
reproduced from [s^]) 
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This quantity, which can in principle be probed by a 
nuclear magnetic resonance (NMR) experiment, was cal- 
culated for the OT approximant. FiglM] shows the six 
different curves obtained by averaging the local suscepti- 
bilities for all sites of a given value of z. As can be seen, 
all of the sites have a low frequency tail, but the different 
types of site have sharp peaks at energies that increase as 
z increases. This provides a precise way of studying local 
environments: by tuning the frequency, one can probe 
different families of sites on the quasicrystal. 
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VII. GEOMETRICALLY DISORDERED 
QUASIPERIODIC QUANTUM 
ANTIFERROMAGNETS 



In this section we consider the effects of geometric dis- 
order due to localized phason flips. By "phason flip" we 
will mean a local change of connectivities as shown in 
Fig |25[ where a site in the interior of moves slightly from 
its original position. This has the effect of modifying 
the bonds locally but leaving the rest of the tiling unaf- 
fected. In addition, the bond lengths and angles remain 
unchanged, as the final result of the phason flip is simply 
a local rearrangement of the tiles. Many quasiperiodic 
structures display this type of phason disorder. The pha- 
sons can be quenched (static) in or they may be excited 
at finite temperature (dynamic). In the present context, 
we consider that the phason disorder is quenched into 
the system (the bonds do not fluctuate between configu- 
rations) , and we are interested in the effect of this type 
of disorder on the antiferromagnetic ground state prop- 
erties. 

The consequences of this type of disorder are, as we 
will see, very different from the consequences of the more 
commonly considered models of disordered magnets, in 
which the couplings or the site occupations are assumed 
to be random. In the case of random couplings f59| . 
the antiferromagnetic ground state is weakened progres- 
sively and ultimately destroyed as the amplitude of the 
fluctuations of bond strengths increases. In the case of 
random site dilution [60| . the antiferromagnetic state is 
weakened and a percolation transition can occur. The 
calculations are carried out on samples with weak pha- 
son disorder. By this we mean that the starting point is a 
perfectly ordered approximant of the quasiperiodic tiling 
(the Penrose or the octagonal tiling) and one randomizes 
these by flipping a certain number of sites so as to create 
a set of disordered samples. The amount of phason dis- 
order introduced in this manner is necessarily bounded, 
as opposed to the random tilings described, for example, 
in 1611. 




FIG. 25: (Color online) A single phason flip showing the orig- 
inal (ro) and final positions {vq, connected with rugged links). 
The original sites and bonds are shown in blue and the new 
site and new bonds in red. (Reproduced from [gj) 



As Figl25] shows, the site undergoing a phason jump 
changes its sublattice. This is taken into account in the 
random selection of flipped sites, since the LSW and 
QMC calculations are performed on samples with an 
equal number of spins on each sublattice, Na — Nb- 
Fig. E51 shows a portion of a perfect (deterministic) tiling 
and a typical example of a disordered tiling after a cer- 
tain numbers of phason flips. For a given total number of 
flips Nph, the degree of disorder A is calculated by calcu- 
lating the average overlap distance between the perfect 
sample and the disordered samples. 




FIG. 26: Samples of (up) a perfect tiling and (down) a ran- 
domized tiling. (Reproduced from [gJ) 



The model considered in 
Heisenberg model for spin- 



64| is the unfrustrated 
described hy H — 
JJ2{ij) Si ' Sj fo'^ spins situated on vertices of approx- 
imants of the Penrose tiling. Nearest-neighbor interac- 
tions are antiferromagnetic J > 0, and act between pairs 
of sites that are linked by an edge. The model is disor- 
dered, with purely geometric disorder, i.e. the couplings 
J acting along bonds remain flxed at a constant value. 
This is to be contrasted with a situation where couplings 
are random, for a given flxed connectivity of sites, as 
in the random exchange antiferromagnet on the square 
lattice [Hi]. 

The results for the sample averaged ground state en- 
ergy of the phason disordered Penrose antiferromagnet 
are shown in Fig |27l The values have been normalized 
with respect to the value obtained for the undisordered 
system. The results obtained from LSWT and QMC 
are in good agreement. This also shows the applicabil- 
ity of the linear spin wave approximations to this ran- 
domized system. The smooth curve shows a fit to the 
form eo = e<jis + (eper - edis)e~°''^, with a = 11.16, 
Sdis = —0.6500 and Cper = —0.6429. The asymptotic 
value of the GS energy Cdis represents the average value 
of the ground state energy of maximally randomized Pen- 
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rose approximants, which lies below the GS energy of 
the perfect system. In other words, overall, the nearest 
neighbor bond energy increases in absolute value. This 
indicates that the introduction of phasons tends to in- 
crease quantum fluctuations in the tilings. This is borne 
out in the calculations of the average staggered magneti- 
zation which should correspondingly decrease when dis- 
order increases. This is indeed seen in Figi28l where the 
average staggered magnetization decreases as a function 
of the disorder strength. Figl29] shows the magnetiza- 
tion profile by an intensity profile, with the perfect tiling 
shown on the left and a moderately disordered tiling on 
the right, for comparison. The colors in the figure on the 
right can be observed to be shifted towards the red end of 
the spectrum, indicating globally smaller values of msi- 
The low energy tail of the IDOS (integrated density 
of states) is quadratic, and can be fitted to give a spin 
wave velocity that increases with A as shown in the 
inset of Fig. [28l This indicates that spin wave propa- 
gation is facilitated by the phason disorder, in analogy 
with the problem of quantum diffusion of electrons in 
the tight binding model in quasiperiodic tilings [651. In 
addition, one notices that the peak corresponding to lo- 
calized states at i? = disappear progressively. These 
states arise, as shown in [55| on closed loops of 3-fold 
sites. They are destroyed when a phason flip occurs on 
the one of the participating sites. 
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FIG. 27: The dependence of the ground state energy on dis- 
order A calculated in LSWT. The smooth curve is a fit to an 
exponential form (see text). The system size is N=4414. The 
inset shows the normalized ground state energy as a function 
of increasing disorder A. Points represent values obtained af- 
ter finite size scaling of LSWT and QMC data. (Reproduced 
from [iil) 



In conclusion, geometrical disorder has important ef- 
fects on the magnetic properties of a quasiperiodic an- 
tiferromagnet. The real space properties of the ground 
state staggered magnetization profile are modified sig- 
nificantly. In particular, the hierarchical structure and 
five-fold symmetries are lost. The overall order parame- 
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FIG. 28: Normalized average staggered magnetization depen- 
dence from the phason disorder A. Points represent values 
obtained after finite size scaling of LSWT and QMC data. 
The inset shows the effective spin wave velocity as a function 
of disorder A for several system sizes. (Reproduced from [6^1 
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FIG. 29: (Color online) Color representation of the local stag- 
gered magnetizations in a perfect tiling (left) and a disordered 
tiling (right) showing sites of magnetizations ranging from 
high (blue) to low (red) values. (Reproduced from [64]) 



ter is smaller, due to the fact that there are more quan- 
tum fluctuations compared to the perfect tiling. It was 
also shown in 55| the effective spin wave velocity for low 
energy modes is slightly higher, showing that disorder 
"helps" propagation of spin waves. Real space properties 
such as the spatial modulation of the order parameter 
are smoothed, and as regards the energy spectrum, the 
fluctuations of the density of states are smoothed. De- 
generate states localized on closed loops disappear with 
increasing disorder. Eigenmodes tend to be more spa- 
tially delocalized compared to the perfect tiling. All these 
effects are peculiar to quasicrystals, for example, the de- 
localization of eigenstates due to phason disorder has also 
been observed in phonon [H^I and tight binding electronic 
models 5l|. These results, obtained for systems with an 
upper bound on the disorder, are not expected to neces- 
sarily hold when one allows for unbounded phason disor- 
der. In the limit of strong disorder, antiferromagnets on 
random tilings in high codimension 63[ are also interest- 
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ing systems for future studies. 



VIII. DISCUSSION AND CONCLUSIONS 

The quantum antiferromagnet for quasiperiodic tilings 
without frustration has been described using different 
theoretical methods. Since the antiferromagnetic cou- 
pling was assumed to be the same for all the pairs of 
spins, the complexity of the Hamiltonian in the tilings 
arises solely from the aperiodic geometry. Two tilings 
were considered in particular: the Penrose rhombus tiling 
and the octagonal tiling, both with and without disor- 
der. We compare the results for the antiferromagnets 
in these tilings with the square lattice antiferromagnet, 
which has, overall, the same number of bonds per spin 
and therefore the same classical ground state energy. 

As in the square lattice, the ground states of the tilings 
have long range magnetic order, of the Neel type. Unlike 
the lattice, the local staggered magnetization nisi is spa- 
tially inhomogeneous, and has a value that depends on 
the type of site i. In the perfect undisordered tilings, the 
values of nisi depend primarily on the number of nearest 
neighbors but also on further near neighbor shells. The 
inflation symmetry of the octagonal tiling allows a calcu- 
lation of the ground state properties by an approximate 
renormalization group. The results agree with the results 
obtained by Quantum Monte Carlo for this tiling. The 
main result of these calculations was to obtain the ground 
state energies of the quasiperiodic tilings and show that 
they are higher in energy for the periodic square lattice, 
which has the same number of spins and bonds. This is 
the first indication that the aperiodicity leads to a sup- 
pression of quantum fluctuations. The staggered mag- 
netizations are in keeping with this conclusion, as the 
average value of m^. is higher in the quasiperiodic tilings 
than on the square lattice. 

To investigate these antiferromagnets in more detail, 
a linear spin wave theory was presented, and solutions 
obtained by numerical diagonalization. The results for 
the staggered magnetization distribution agree with that 
found by other methods. The dependence of nig on the 
local environments was investigated in some detail, and 
explained using a Heisenberg star cluster argument. The 
staggered magnetization in real space was seen to be a 
complex but highly symmetric function which is simpli- 
fied when represented in perpendicular space coordinates. 
The Fourier transform of the magnetic state can be de- 
scribed by a quasiperiodic peak indexing and is related 
to the original nuclear structure factor by a generalized 
(i.e. in dimensions greater than the physical d = 2) wave 
vector shift. This feature, characteristic of standard two 
sublattice antiferromagnets, has its counterpart in real 
three dimensional magnetic quasicrystals, as measured 
in neutron scattering experiments. 

The spin excitations on the tilings have complex spa- 



tial structure, and show multifractal scaling. Their spa- 
tial extension is highly dependent on the energy, with the 
states of low energy being the most extended as seen by 
the scaling of the IPR (Inverse Participation Ratio) with 
system size. The correlations are long ranged. Spin wave 
velocities, as found from fitting to the low energy end of 
the magnon spectrum decrease with increasing complex- 
ity of the tilings. The local dynamical susceptibility is 
shown to have sharp peaks as a function of frequency, 
and depends sensitively on the local environment. Thus, 
magnetic response could be a probe for self-similarity of 
quasiperiodic structures. The entanglement properties of 
such systems should present interesting differences with 
respect to the periodic lattice [65|] . Interesting nontrivial 
behavior of entanglement entropy has been found in one 
dimensional aperiodic spin models (66| but two dimen- 
sional quasiperiodic spin models have not been studied 
so far. 

Finally, the effects of phason flip disorder were de- 
scribed. The ground state increases as disorder increases, 
in the weak disorder limit. This leads to the conclusion 
that phason disorder makes the system more classical ~ a 
conclusion also borne out by the fact that magnetizations 
globally increase in value. This is the opposite to the be- 
havior usually found in crystals, where disorder tends to 
diminish the antiferromagnetic order parameter. Wave- 
functions tend to become more extended, when phason 
disorder is increased - this behavior was also observed in 
the tight-binding model for electrons in quasicrystals. 

The two tilings studied are the best known of the two 
dimensional quasicrystals. Other tilings that could be 
investigated, in order to better understand the quantum 
fluctuations in aperiodic systems and the effects of dis- 
order, are the 2d Rauzy tilings, and the rhombus tilings 
of Vidal et al. We have discussed unfrustrated models, 
on bipartite tilings. However in real quasiperiodic com- 
pounds, it is more likely that the structures are not bipar- 
tite, leading to frustrated antiferromagnetic interactions. 
This type of problem was recently discussed for classi- 
cal spins 67[ but has not been investigated for quantum 
models. 

In conclusion, quasiperiodic antiferromagnets are a ex- 
ample of highly complex order in a magnetic system. Ex- 
perimentally, the Heisenberg model may be possible to 
realize in thin films grown by deposition on quasicrystal 
surfaces. Cold atoms provide an alternative laboratory 
to test the theoretical predictions. 
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